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Abstract 

The heavy ion probability for continuum e + e~ pair production has been calculated to all orders 
in Za as a function of impact parameter. The formula resulting from an exact solution of the 
semiclassical Dirac equation in the ultrarelativistic limit is evaluated numerically. In a calculation of 
7 = 100 colliding Au ions the probability of e + e~ pair production is reduced from the perturbation 
theory result throughout the impact parameter range. 
PACS: 25.75.-q, 34.90. +q 
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I. INTRODUCTION 



The impact parameter dependence of heavy ion e + e~ pair production to all orders in 
Za is of direct relevance to experiments carried out at the BNL Relativistic Heavy Ion 
Collider (RHIC) and possible at the CERN Large Hadron collider (LHC). In the usual 
method of scattering theory, the impact parameter is formally integrated over to turn the 
momentum transfer of the amplitude into that of the cross section via a delta function[l|. 
However a standard method of triggering measured events at RHIC and LHC is with zero 
degree calorimeters, devices that respond to neutrons from Coulomb dissociation of the 
beam ionsj^j]. The physical process of a measured pair production event then also includes 
simultaneously the Coulomb dissociation of the ion nucleus. Assuming independence of the 
two processes (pair production and Coulomb dissociation) the probability of an event is 
then the product of the probability of pair production times the probability of Coulomb 
dissociation. Predicted cross sections are then appropriately constructed from the integral 
of the product of the pair production probability times the dissociation probability. I will 
confine myself in this paper to evaluation of the impact parameter dependent probability 
for heavy ion production of e + e~ pairs. Combining the probability of the particular e + e~ 
phase space acceptance of physical detectors with probability of Coulomb dissociation will 
be left for future treatment. 

The impact parameter dependent calculation of heavy ion e + e~ pair production to all 
orders in Za described in this paper is based on techniques previously developed. That de- 
velopment began with the realization that in an appropriate gauge[3(, the electromagnetic 
field of a relativistic heavy ion is to a very good approximation a delta function in the direc- 
tion of motion of the heavy ion times the two dimensional solution of Maxwell's equations 



in the transverse direction 



4| , and that an exact solution of the appropriate Dirac equation 



could be obtained for excitation of bound-electron positron pairs exhi 
perturbation theory of a little less than 10% for Au + Au at RHIC 



Ditine a reduction from 

3. 



An analytical solution of the Dirac equation was subsequently obtained 0, Q, Q] for the 
analogous case of continuum e + e~ pair production induced by the corresponding counter- 
moving delta function potentials produced by heavy ions in an ultrarelativistic collider such 



as RHIC. It was noted at the time 
perturbation theory result. 



3, |8| that the exact solution appeared to agree with the 



B, y 



Several authors then argued [9|, [10|, llll] that a correct regularization of the exact Dirac 



equation amplitude should lead to a reduction of the total cross section for pair production 
from perturbation theory, the so called Coulomb corrections. The first analysis was done in 
a Weizsacker- Williams approximation 9]. Subsequently, Lee and Milstein computed jiol. 
the total cross section for e + e~ pairs using approximations to the exact amplitude that 
led to a higher order correction to the well known Landau-Lifshitz expression jl^]. One 
might also have noted that it is well known that photoproduction of e + e~ pairs on a heavy 
target shows a negative (Coulomb) correction proportional to Z 2 that is well described by 
the Bethe-Maximon theory {13]]. In fact some years ago Bertulanni and Baur]14] applied 
the structure of the Bethe-Maximon theory to obtain an impact dependent probability for 
heavy ion reactions with a negative Coulomb correction. I will comment on their results 
and compare them with the results of the present work in Section V. 



In a previous paper I have tried to exp 



icate the Lee and Milstein approximate results 



and argued their qualitative correctness [151]. Subsequently I undertook the full numerical 



calculation of electromagnetically induced ultrarelativistic heavy ion electron-positron pair 



production 16]. That evaluation of the "exact" semi-classical total cross section for e + e 
production with gold or lead ions showed reductions from perturbation theory of 28% for the 
CERN Super Proton Synchrotron (SPS) case, 17% for RHIC, and 11% for LHC. For large 
Z no final momentum region was found in which there was no reduction or an insignificant 
reduction of the exact cross section from the perturbative cross section. 

In the present paper I reformulate the approach used for my previous cross section cal- 
culations in order to calculate corresponding impact parameter dependent probabilities. 

II. IMPACT PARAMETER DEPENDENT PROBABILITIES 

For production of continuum pairs in an ultrarelativistic heavy ion reaction one may 
work in a frame of two countermoving heavy ions with the same relativistic 7, and the 
electromagnetic interaction arising from them goes to the limit of two 5 function potentials 

V(p, z, t) = 8{z - t)(l - a z )A-(p) + 5(z + t)(l + a z )A + (p) (1) 

where 
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I calculate the number weighted probability (or number operator) < N(b) > for producing 
e + e~ pairs at some impact parameter b 

< N(b) >=±nP n (b) = / "ff l d3q \M(p,q)\* (3) 

n=l J ( 27r ) % e <? 

with the previously derived exact semiclassical amplitude for electron-positron pair 
production (], 7, 8] written in the notation of Lee and Milstein[lo| 

M(p,q) = J ^exp[zk-b]A<(k)F B (k)F A (q ± + p ± -k), (4) 

where p and q are the four-momenta of the produced electron and positron respectively, 
P± = Po ± Pz, q± = Qo ± Pzi 7± = 7o ± 7z, ot, = 707, k is an intermediate transverse 
momentum transfer from the ion to be integrated over, 

km\ \ -f \ a ' (k- p±) +7o^ , . 

M{k) = u{p) — — 7_M (-q) 

—P+Q- — (k — pj_)^ — m 2 + le 

. -/ \ -a • (k - q±) + 7 m 

~P-Q+ ~ (k — q±) — + 



and the effect of the potential 



jq. (1-2) is contained in integrals, F B and Fa, over the 



transverse spatial coordinates 0, 3,0], 

F(k) = J d 2 p exp[-ik- p]{exp[-iZalnp] - 1}. (6) 

F(k) has to be regularized or cut off at large p. How it is regularized is the key to under- 
standing Coulomb corrections as I will review in Section III. 

In Section III I will show how a properly regularized transverse potential in Eq. (6) 
necessitates a numerical integration over p to obtain F(k). With the resultant F(k) the 
integration over k in Eq. (4) is then not amenable to the usual Feynman integral techniques 
and must be carried out numerically. For numerical convenience I have chosen to express k 
in Cartesian co-ordinates and rewrite Eq. (5) 

a x k x + a y k y - a p - 7 m 

M(k) = u{p) v -f - — - 7 _u -g 

—P+1- — (k — P_l) ~~ m + te 
, _/ \ -a x k x - a y k y + a • q ± + 7o™ , v 
~P-Q+ ~ (k — q±J — w + ^e 

The expression for the amplitude M(p,q) Eq.(4) then becomes 

M(p,q) = u(p)[I px a x + I py a y + (-a • p± + 7 m) J p ]7_w(-g) 

+M(p)[-/^ax - ^j/tt?; + (a ■ q± + 7 m) J g ]7 + n(-g), (8) 



where, letting b define the x-axis, 

/ -— [ex»\ik b\dk / ^( k )^(q± + P±-k) k x dk y 
Ipx ~ (2vr) 2 J eXV[t kx 6J dkx J - V + Q _ — fk — d i ) 2 — m 2 

i» = ^ I «p[< *■ & ] / t: v do) 

V = (2^2 / ex P^ *■ & ] / -^-^ t ; ~ (12) 



- (k - p± 


) 2 — m 2 


Ffl(k) F A (q ± + P± - 
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) 2 — m 2 
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|2 _ m 2 
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2 — m 2 



J P = ~MjA% I exp ^ kx & ] ^ / — ~ 7V, — ^2 — z^r ( 13 ) 



(2k) 

After squaring, _g over sp.n states, and taking traces with the aid of the computer 
program FORM 17 I obtain the expression for the amplitude squared 

\M(p,q)\ 2 = P+ q4(m 2 +pl)\J p \ 2 + \I px \ 2 + \I py \ 2 - 2p x Re(J p P px ) - 2p y Re(J p P py )} 
+ p^q + [(m 2 + q\)\J q \ 2 + \I qx \ 2 + \I qy \ 2 - 2q x Re(J q P qx ) - 2q y Re( J q I*J] 
+ 2[(m 2 +p 2 ± )(q x Re(J p r qx ) + q y Re(J p P qy )) 

+ (m 2 + q 2 ± )(p x Re(J q P px ) + p y Re(J q P py )) 

+(p± ■ q± - m 2 )(Re(I px r qx ) + Re(I py P qy )) 

-(m 2 +pl)(m 2 + q 2 ± )Re(J p J* q ) 

-(PxQy + p y q x ){Re(I px Ily) + Re(IpyI* x )) 

-2p x q x Re(I px I* x ) - 2p y q y Re(I P yI* y )}. (15) 



There is an apparent numerical difficulty in evaluating Equations (9)-(14) due to the 
oscillating factor exp[ik x b]. Of course in the b = limit this factor is absent, and so we 
will first investigate this numerically more tractable case. The general case of non-zero b 
will then be addressed by a technique involving piecewise analytical integration. 

The derived impact parameter dependent amplitude squared, \M(p, q)\ 2 , is not simply 



the square of the amplitude for the excitation of a particular (corre 



pair 



18(. However, as has been discussed in the literature 



20 



ated) electron 



21 



22 



23|, 



24 



positron 



25|, 



26] 
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the number operator for a particular (uncorrelated) electron or positron can be constructed 
by integrating over the positron or electron momenta respectively, and likewise the number 
operator for total pair production. In previous articles Q I have also discussed these 
matters in more detail. 

For the present I will limit myself to calculating the number operator < N(b) > for total 
pair production, and from it the total pair production cross section <jt 

< N ®>=iw0r M[2 (i6) 

Note that <jt corresponds to a peculiar type of inclusive cross section which we should call 
the "number weighted total cross section" , 

a T = d 2 b < N(b) >= / d 2 b ]T nP n (b), (17) 

n=l 

in contrast to the usual definition of an inclusive total cross section 07 for pair production, 

„ 00 

ai = J d 2 bY, p niP). (18) 



n=l 



III. EIKONAL, EXACT, AND PERTURB ATIVE CASES 

If one merely regularizes the transverse integral Eq. (6) at large p one obtains[?j, 0] apart 
from a trivial phase 



F M = T^z- (19) 



AwaZ 
~k 2 

Then all the higher order Za effects in M(p, q) are contained only in the phase of the 
denominator of Eq. (19). It was noted that as Z — > this photon propagator leads to the 
known perturbative result for M(p, q) if F(k) is modified to 

b W = (P + ^/7 2 ) 1 -^- 1 } 

In this approach a lower k cutoff at some u/j has to be put in by hand to obtain dependence 
on the beam energy and to agree with the known perturbative result in that limit. 

Rather than put in the uj 2 /'-/ 2 cutoff by hand, my physically motivated ansatz is to apply 
a spatial cutoff to the transverse potential x(p) (which has been up to now set to 2Za In p) 



in order to obtain an expression consistent with the perturbation theory formula 27|, [28] 



in the ultrarelativistic limit. In the Weizsacker- Williams or equivalent photon treatment 



of electromagnetic interactions the effect of the potential is cut off at impact parameter 
b ~ where 7 is the relativistic boost of the ion producing the photon and ui is the 

energy of the photon. If 

X (p) = r dzV(y/z* + ?) (21) 
and V(r) is cut off in such a physically motivated way, then 11] 



where 



V(r) = - Za ^M-^A,B/l] (22) 



P+ +Q+ P-+ q- /nn\ 

= — 7, — ; u B = — 2 — ( 23 ) 



with uo A the energy of the photon from ion A moving in the positive z direction and ujb the 
energy of the photon from ion B moving in the negative z direction. Note that we work 
in a different gauge than that used to obtain the original perturbation theory formula, and 
thus our potential picture is somewhat different. The transverse potential will be smoothly 
cut off at a distance where the the longitudinal potential delta function approximation is no 
longer valid. 

The integral Eq. (21) can be carried out to obtain 

X (p) = -IZuK^puiA^h), (24) 

and Eq. (6) is modified to 

F A>B (k) = 2tc J dppJ (kp){eyLp[2iZ AjB aK (puj AiB /'j)] - 1}. (25) 

F A (k) and F B (k) are functions of virtual photon u A and uj b respectively. The modified 
Bessel function 

K {pujh) = - ln(p) - ln(aV2 7 ) - 7e (26) 

for p <ti (7 e is the Euler constant), and cuts off exponentially at p ~ 7/0;. This is the 
physical cutoff to the transverse potential. One may define £ = kp and rewrite Eq.(25) in 
terms of a normalized integral I AjB {^k/uj) 



Fa,b(V = A -^^I A , B (lkM (27) 



where 

IaM^M = ttt^ / dZZJ (Z){exp[2iZ AtB aK (fr/>yk)) - 1}. (28) 




FIG. 1: Top: Relative magnitude squared of virtual photon source; solid line exact; long dashed 
line eikonal or perturbative (see text). Bottom: Corresponding phase of virtual photon source; the 
perturbative phase (not shown) does not vary with k^/u. 

Note that Fa,b is equal to ^tkZa.bI^ 2 times a function of (jk/u). The form of the integral, 
Eq. (6) without the cutoff (here with added arbitrary phase constants consistent with 

nn 

Ko(^uj/jk) for small £, Eq. (26)) was solved in closed form[7|, |8J 

W 7 * /w) = r(iaz) (i + ^AV) 1 -^- (29) 



I will refer to this form as the eikonal. In the limit as Z — > of I^^^k/uj) goes to 



—i 



lA, B {lk/uj) = 2 , (30) 

1 + or/ £r 7 Z 

leading to the familiar perturbation theory form 

F ^ (k) " k* + uW (31) 

The top panel of Figure 1 displays the results of numerical calculation of \I{k^f/uj)\ 2 for 
Z = 79 and in the perturbative limit. The curves that go to zero on the left are \I{k^(/uj)\ 2 . 
The curves that go to zero on the right are \I{k^f/uj)\ 2 multiplied by (u/kj) 4 to exhibit the 
reduction at high Z on the left of the curve. {u / k^) A \I {k~f / u)\ 2 goes to one as kj/uj goes to 
zero in the perturbative case; it goes to a reduced constant value as kj/u goes to zero for 
Z = 79. The bottom panel exhibits the phase of the full expression F(k). The perturbative 
phase (not shown) is constant. Both exact and eikonal phases are identical at high k^/uo 
and diverge by 22 degrees for this Z = 79 case as they approach constants at low k^/u). 

IV. THE b = LIMIT 

Since the present calculations, based on the ultrarelativistic limit, employ a somewhat 
different approach than previous perturbative calculations, but should agree with them 
in the perturbative limit, I have deemed it prudent to do some comparison of numerical 
results. As a first step I computed the b = limit, both as a test of the method and general 
structure of the computer code as well as for physical interest. Figure 2 shows a pertubation 
theory calculation of the energy spectrum of produced positrons at b = compared with a 



previous perturbation theory calculation of Hencken, Trautmann, and Baur 30]. The same 
dipole form factor has been utilized for the ions with A = 83 MeV, and the curves have 
been divided by (Za) A . The heavy ions are from Au + Au at RHIC, and the energy of the 
electrons has been integrated over. Except at the lowest energies the agreement is good, 
and the integrated probability of the present calculation P°(0) = 1.64 is in good agreement 
with the previous calculation P°(0) = 1.6. The integrated numbers have not been divided 
by (Za) 4 as the displayed curves were. 

Figure 3 shows the effect of higher order contributions. My perturbation result of Figure 
2 is repeated for reference. The solid line is the result of the exact calculation. It is clearly 




FIG. 2: Total probability as a function of energy of the positron. The curves have been divided 
by (Za) 4 . Short dashes: present perturbation theory; Dotted line: previous perturbation theory 



calculation 



30] 



reduced from perturbation theory for the lowest to highest values of the energy range; the 
exact energy integrated P(0) = .94 = .57P°(0). Only slightly lar ger in magnitude is the 
result of an eikonal calculation where the energy integrated P E (0) = 1.03 = .63P°(0). 

The overall reduction from perturbation theory is almost as large in the eikonal case as it 
is in the exact calculation, and the two curves in Figure 3 do not diverge much at any energy 
value. This would suggest that similar effective physics is involved in the two calculations. 
Recall from Figure 1 the reduction in magnitude of the virtual photon source for the exact 
calculation as compared to the eikonal (which is identical in magnitude to perturbation 
theory) for k^/uj < 50. For k > 100 the magnitude of P(k) goes over into the eikonal 
or perturbative result, while the phase of F(k) is identical to the eikonal. It is in this large 

10 



0.01 




Positron Energy E (MeV) 



FIG. 3: Probabilities as in Fig.2. Short dashes: present perturbation theory; Solid line: exact); 
Long dashes: eikonal 

k region, where F(k) becomes identical the eikonal in both magnitude and phase, that most 
of the contribution for b = comes. Test calculations show that 80% of the contribution 
to the total probability at b = comes from k > 100 and that less than a tenth of a 
percent comes from k < 10 co/j. It is not the region of reduced magnitude that dominates 
for the exact case b = 0, but the region of rotating phase, and that rotating phase is what 
reduces the cross section from perturbation theory. The eikonal F(k) is a fair approximation 
to the exact at b = 0, but as we will see in Section V this is not true for larger b. 

It is worth emphasizing that calculations labeled perturbative, exact, and eikonal differ 
only in the expressions used for Fa,bQ*). The analytical expression Eq. (31) is used for 
perturbative calculations, and Eq. (27) and Eq. (29) for the eikonal calculations. The exact 
calculations makes use the expression Eq. (27) and Eq. (28), which must be evaluated 
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numerically, but only once for each Za,b of interest. 



V. b DEPENDENT PROBABILITIES AND THE TOTAL e+e CROSS SECTION: 
Au + Au AT RHIC 

With the addition of impact parameter dependence, we are faced with a nine dimensional 
integral for the total cross section as compared to the seven dimensional integral in the previ- 
ous section, or the seven dimensional integral in the representation of cross sections without 
impact parameter [3]. Although the usual method of evaluation, e.g. in perturbation theory, 
is via Monte Carlo, I have chosen to do the multi dimensional integral directly on meshes 
uniform on a logarithmic scale in each radial momentum dimension and on a logarithmic 
scale in impact parameter. It was computationally tedious but possible to carry out the 
calculations to some rough but significant accuracy without using Monte Carlo because the 
integrands are very smooth and smoothly goes to zero at both high end and low end of the 
momentum ranges and impact parameter range. 

Apart from the fact that impact parameter dependence adds two more dimensions to 
calculations (the magnitude of b and the angle of the outgoing pair relative to b), evaluation 
of the square of the pair production amplitude, Eq. (15), is considerably more difficult for 
non-zero impact parameter b due to the factor exp[i k x b] in Eq. (9-14). For increasing values 
of b the integral over k x oscillates rapidly and presents a challenge to numerical integration. 
Furthermore, while the oscillation is on a linear scale in k x , experience has shown that the 
natural scale for integration over momenta in this electromagnetic problem is logarithmic; 
we want to choose numerical mesh points on a logarithmic scale. I have chosen to work in 
Cartesian coordinates because after a numerical integration over the comparative smooth 
variation in k y I am able to make use of a relatively simple piecewise analytical method to 
integrate over k x . I do the integration over k x in the manner described in Appendix A. 

Calculations were again carried out for the RHIC case of Au + Au at 7 = 100. For 
the perturbation theory and eikonal probabilities I was able to do a rough check of my 



29J . Figure 4 shows 



computer code by comparing with another set of published calculations 
the comparison of cross section contributions derived from the present calculations with 
those derived from the b dependent probabilities of Hencken, Trautmann, and Baur. The 
logarithmic derivative was chosen so that contributions to the integrated cross section go as 
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FIG. 4: Impact parameter dependence of contribution to total cross section. Dashed line, perturba- 
tion theory; long dashed line eikonal. Comparable results are shown derived from the calculations 
of Hencken, Trautmann, and Baur 29j]: Dotted line, perturbation theory; Dot dashed line, eikonal. 



the area under the curves on this natural log scale plot. No form factor was utilized in either 
calculation as its effect should be small for total cross section probabilities at these impact 
parameters. Agreement is reasonably good, considering the large grids used in the present 
calculations. Note that that for both sets of calculations the eikonal solution is less that 
perturbation theory for low impact parameters, and then crosses over to be greater than 
perturbation theory at about b = 3000 fm and larger. For small impact parameters there 
is fair agreement, but there is an incipient divergence between the curves for the present 
calculations and those of Hencken, Trautmann, and Baur where the latter cut off. However, 
one of the authors had warned me not to trust their results far out in b too much because 



the Fourier transform used was not as good there 3l| 
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Impact Parameter b (fm) 



FIG. 5: Impact parameter dependence of contribution to total cross section. Dashed line, pertur- 
bation theory; long dashed line eikonal; solid line, exact. 

In Figure 5 one can observe the comparison of the exact calculation with perturbation 
theory and the eikonal. For all impact parameters the probability of pair production is 
smaller in the exact calculation than it is in perturbation theory. The small impact parameter 
agreement between the exact calculation and the eikonal extends the mechanism discussed 
in the b = case, the reduction from pertubation theory due to the similar rotating phase of 
F{k) for dominant large values of k seen in Fig. 1. For larger impact parameters providing 
the dominant contribution to the cross section, a significant contribution comes in the region 
k < 10 toj^j where the magnitude of F(k) is reduced from perturbation theory and the 
eikonal. This demonstrates how a correct physical cutoff process preserved the Coulomb 
corrections at high b and thus leads to Coulomb corrections in the total cross section. 

The calculations can be integrated over impact parameter to compare with my previous 
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FIG. 6: Impact parameter dependence of contribution to total cross section. Present calculation: 
Dashed line, perturbation theory; Solid line, exact. Comparable results are shown derived from 



the calculations of Bertulani and Baur[l4l]: Dotted line, perturbation theory; Long dashed line, 
Coulomb corrected. Comparable perturbative results derived from the formula of Lee, Milstein, 
and Serbo 32| are represented by the dot-dashed line. 



calculations fig] done not in the impact parameter representation. That the present inte- 
grated perturbation theory calculation of 34.6 kbn (kilobarns) is in such good agreement 
with my previous calculation of 34.6 kbn must be considered fortuitous, given the relatively 
crude nature of the impact parameter calculation. The present integrated exact calcula- 
tion of 29.4 kbn is in fair agreement with my previous 28.6 kbn. The integrated eikonal 
calculation of 35.5 kbn should in principle be identical to the 34.6 kbn perturbation the- 
ory calculation, but like the present exact calculation is slightly larger than my previous 
published non impact parameter calculation. 
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At this point it is interesting to also compare the present results with two previous ap- 
proximate analytical expressions for impact parameter dependence. In their classic Physics 
Reports article, Bertulani and Baur 14| presented an impact parameter formula that included 
a Coulomb correction. Figure 6 shows a comparison of the cross section contribution derived 
from their probability formula with the present results. Bertulani and Baur integrated over 
their stated range of validity, from l/m e to ~ 7 2 /m e , to obtain a Landau-Lifshitz perturba- 
tive result with a Coulomb correction. For symmetric heavy ions, the Coulomb correction 
that they obtained was of the same form as the result later obtained by Serbo et al. 9| and 
by Lee and Milsteinfiol. \v\ . However, since Bertulani and Baur had postulated averaging 
the result of a Coulomb correction for the target with one for the projectile, they obtained 
half of the more correct later result based on adding the lowest order Coulomb correction 
for both target and projectile. The overall shapes of the Bertulani and Baur results differ 
from the present results, but their Coulomb correction reduces the probability throughout 
the impact parameter range. 

Shortly after obtaining their above mentioned Coulomb correction results for the total 
cross section t], Q, 11] Lee, Milstein,and Serbo presented an analytical form for the the 
impact parameter dependence}^] that differed substantially from that of Bertulani and 
Baur, especially in the intermediate range from l/m e to r y/m e . Figure 6 also displays 
the impact parameter contribution to the cross section obtained from their pertubative 
expression for the probability. The shape of the curve seems qualitatively more in agreement 
with present results. As in the case of Bertulani and Baur, the treatment is considered valid 
for values of impact parameter larger than l/m e (386 fm), and in both cases integration 
of the contribution from their expression is in agreement with the leading ln 3 (7 2 ) term of 
Landau and Lifshitz 121] . Lee, Milstein, and Serbo did not present a specific form for the 
impact parameter dependence of the Coulomb correction, but as noted above, their previous 
work indicates a negative Coulomb correction on average twice that of Bertulani and Baur 
over the indicated parameter range. 



VI. SUMMARY AND DISCUSSION 

Calculated exact total probabilities for heavy ion e + e~ pair production exhibit a reduction 
from the probabilities calculated in perturbation theory throughout the full range of impact 
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parameters. 

In principle the impact parameter approach to calculating exact e + e~ pair production 
probabilities is suited for combining with an impact parameter dependent mutual Coulomb 
dissociation calculation [2] to be able to com par e with e.g. data such as was obtained with the 
Solenoidal Tracker at RHIC (STAR) setup [33]. In this paper I have presented calculations 
only for the total e + e~ pair production probabilities. Still to be done is a sufficiently accurate 
calculation of the high transverse momentum slice of data seen by STAR to be combined with 
a Coulomb dissociation calculation for the zero degree calorimeter acceptance. Also since the 
present approach is strictly speaking valid only when either the positrons or electons have 
been integrated over, and in the STAR case both electron and positron are constrained to be 
in the high momentum slice, the present approach to e + e _ production is not exactly valid for 
the STAR case. At present the best one can do is observe that the present method is valid for 
both uncorrelated positrons and electrons of all momenta, and ignore the correlations. The 
effect of correlations averages to zero, but some estimate of individual magnitudes would be 
useful. 
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APPENDIX A: PIECEWISE ANALYTICAL FOURIER INTEGRATION 

Consider the integral over k x (I will drop the subscript x for now). It contains a smoothly 
varying part which I will call f(k) times the rapidly varying coefficient exp[i kb] 

f(k)exp[ikb}dk. (Al) 

One can use integration by parts to transform I to 

[l-exp[ikb])<tf(k) 
ib dk 
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dk, (A2) 



where the term 1 is inserted for numerical convenience. One now rewrites J in a form 
suggesting piecewise analytical integration from mesh point to mesh point 

fc, max r k i+1 (l - exp[i df(k) 



I = 



^ Jki ib dk { ' 

Taking the lowest order (constant) approximation to the derivative over the interval one has 

j = k y( f(h + i)-f(h) \ /-fci+i (l-exp[ikb]) dk 
^ V k i+ i -h J Jk r ib 



H+l 



k ^/ f(k i+1 )-f(h) \/ k i+1 -h exp[ik l + 1 b]-exp[ik i b] \ 

fc M h+i-ki )\ ib v )■ 1 } 

This is the lowest order expression for piecewise analytical integration. 

One can obtain the piecewise analytical expression good to the next leading order by 
taking three point Lagrange interpolation for f(k) between ki and k i+i and then taking the 
derivative. One obtains 

^imax 1 / Q 



"-imax 1 / LJ 

1 = E - h (Mh+i-h) + -(k\^-k\) 



- exp[. k l+1 b] ( ^± + -)+ exp[* ki b] (— ^ + ) , (AS) 

where 

A _ f f(h) (fcj+i + ^+2) | f(k i+ i) (kj + k i+2 ) | f(k i+2 ) (ki + fcj+i) \ , Ag s 
V(/ci - fcj+i) (fcj - fc i+2 ) (fcj+i - fcj) (k i+ i - k i+2 ) (k i+2 - h) (k i+2 - ' 



and 



f(h) f(k i+ i) f (ki+2 1 



V (ki - k i+1 ) (ki - k i+2 ) (k i+1 - ki) (k i+1 - k i+2 ) (k i+2 - ki) (k i+2 - k i+1 ) J ' 

The expression Eq. (A5-A7) exhibits numerical difficulties for small impact parameters. 
Therefore, for impact parameters in the range below the electron Compton wavelength 
(about 386 fm.) I have utilized Eq. (A4), and I have utilized Eq. (A5-A7) for higher impact 
parameter values. 
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